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ABSTRACT 

A general prescription to solve the Dirac equation in the presence of CP- 
violating electroweak bubble wall is presented. The profile of the bubble wall 
is not specified except that the wall height is mo and zero deep in the broken- and 
the symmetric-phase regions, respectively, where mo is a fermion mass given by 
the Higgs-vacuum-expectation value and the Yukawa coupling. The CP-violating 
effects are evaluated by regarding CP-violating part of the bubble wall as a per- 
turbation to CP-conserving solutions. The basic quantity, Rr-^l — Rr^l, which 
would contribute to the cosmological baryon asymmetry, is estimated for some 
typical profiles of the wall, where Rr^l{Rr^l) is the reflection coefficient of 
right-handed chiral fermion (anti- fermion). 
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1. Introduction 



It is one of the most challenging problems in particle physics and cosmology 
to explain the baryon asymmetry of the universe. In order to generate the baryon 
asymmetry, the basic theory should at least satisfy the famous three conditions'^' 
: there exist (i) baryon-number violation, (ii) CP nonconservation and (iii) out-of- 
equilibrium processes. It has been recognized that the baryon number is violated 
in the standard model through change of gauge- field configurations due to the axial 

[^] 

anomaly . At very high temperature, the baryon number is badly broken even to 
erase any primordial asymmetry, while the transition rate becomes negligible at 
low temperature. Thus there is a potential possibility in the standard model to 
produce the baryon asymmetry if the electroweak phase transition at a temperature 
T ~ 100 GeV is of first order, which guarantees the above condition (iii), and if CP- 
nonconserving interactions are effective. It is known that the phase transition in 
the standard model may be weakly of first order for light Higgs boson. However, the 
experimental lower bound on the Higgs-boson mass (~ 60 GeV) seems to contradict 
the requirement to forbid anomalous processes induced by the sphaleron after the 

[3] 

electroweak phase transition . Furthermore it seems that the CP violation in 
the Kobayashi-Maskawa scheme'^' is too feeble to generate the observed baryon 
asymmetry. 

A simple generalization of the standard model would be to introduce another 
Higgs doublet . In two-Higgs-doublet models, including the supersymmetric ex- 
tension of the standard model, there could appear a new CP- violating phase in the 
Higgs sector. It is expected that the new CP violation could make baryon asymme- 
try large enough to be comparable with the observed value. Moreover it is argued 
that the constraint on the mass of the lightest Higgs boson to avoid washing out 
of the baryon asymmetry after the electroweak phase transition is not so stringent 
that the Higgs boson could be as heavy as about lOOGeV . It is quite possible 
that the phase transition may be of first order because of many degrees of freedom 
in the Higgs sector'^'. In a previous paper, three of the authors ( K. F., A. K. and 
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K. T. ) have shown that the electroweak phase transition is clearly of first order in 
a two-Higgs-doublet model with a constraint that the Higgs potential contains no 
quadratic terms in the symmetric phase'*''. 

In this paper, we study fermion propagation in CP- violating bubble-wall back- 
ground, which would play an important role to generate the cosmological baryon 

[7] 

asymmetry . We give a general prescription to solve the Dirac equation in the 
presence of CP- violating bubble wall* For simplicity, we take a single species of 
fermion, e.g., the top quark. If the electroweak phase transition is actually of first 
order, supercooling will occur and bubbles of the broken phase will be generated 
in false vacuum. The bubbles will expand with some velocity and the space will be 
filled with the true vacuum. Since the bubble scale is macroscopic, the profile func- 
tion of the bubble wall may be assumed to depend on a single spatial coordinate 
in the rest frame of the bubble wall. Far from the wall, deep in the broken-phase 
region, the Higgs field develops a constant vacuum expectation value. We assign 
a complex mass to the fermion as a function of m{z) — mji{z) -\- imj{z), where 
z is the coordinate perpendicular to the wall. The real part ttir^z) asymptotically 
behaves such that m^{z) as z ^ — oo ( symmetric phase ) and m^{z) mg 
as z +00 ( broken phase ), where mo is the fermion mass. The imaginary part 
mj{z) gives rise to CP violation if mi{z)/mR{z) is not a constant. 

We do not specify any explicit form of the profile function 777.(2;) and obtain 
general expressions for the refiection and the transmission coefficients. The CP- 
violating effects are evaluated by regarding the CP- violating term as a perturbation 
to CP-conserving solutions. Namely, we first solve the CP-conserving Dirac equa- 
tion. Subsequently we solve the full equation to the first order of the CP-violating 
term by means of the Green's function''^'. Our method corresponds to the distorted 
wave Born approximation ( DWBA ). The reflection coefficient R and the transmis- 
sion coefficient T for the fermion incident from the symmetric phase region to the 
wall satisfy the unitarity R + T — 1 to the first order of perturbation theory. The 

7k- For a special bubble- wall profile of the kink type with no CP violation, the Dirac equation 
has been studied by Ayala et all^^ . 
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basic quantity, AR = Rr-,l — Rr-^l, which would be important to explain the 
cosmological baryon asymmetry , is estimated for some typical profiles of the wall, 
where Rr^l ( Rr-^L ) is the reflection coefficient for right-handed chiral fermion 
(anti-fermion). The absolute value of Ai? is found to decrease rapidly with the in- 
crease of E* = E'^ — py, where E and pr are energy and transverse momentum 
of the incident fermion, respectively. It seems that |Ait!| decreases more rapidly 
with the increase of E* for thicker walls. We expect that our prescription would 
serve to construct new models to generate the cosmological baryon asymmetry in 
the electroweak phase transition of first order. 

In the next section, we will present a general prescription to solve the Dirac 
equation in the background of CP- violating bubble wall. The reflection and the 
transmission coefficients are obtained from the solutions. In Section 3, the results 
in Section 2 are apphed to the special case where the real part of the bubble-wall 

rgj 

profile is of the kink type. In this case we estimate Ai? numerically by assuming 
typical forms for the imaginary part of the bubble-wall profile. The final section is 
devoted to conclusions. Details of calculations will be given in Appendix. 

2. DWBA to CP- Violating Dirac Equation 

2.1. Dirac equation and Ansatz 

Here we solve the Dirac equation in the bubble-wall background with CP vio- 
lation. For simplicity, we consider one-flavor model described by the lagrangian, 

C = i^LlHh + ^RiHr + {fi^Li^R(t> + h.c). (2.1) 

In the vacuum, near the first-order phase transition point, (0) may be x-dependent 
field, so that we put 

m{x) = -fmx), (2.2) 
where m{x) is in general complex-valued and we neglect the time dependence. 
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If the phase of m{x) has no spatial dependence, it is removed by a constant bi- 
unitary transformation. We are interested in x-dependent phase which cannot be 
transformed away. The Dirac equation to be solved is 

i0ip{t, x) - m{x)PRip{t, x) - m*{x)PLip{t, x) = 0, (2.3) 

where Pr = (1 +75)/2 and Pl = (1 — 75)72. In the bubble- wall background whose 
radius is sufficiently large, m{x) can be regarded as a function of only one spatial 
coordinate, so that we put m{x) — m{z). 

To solve (2.3), we take the following Ansatz: 

i^it, x) = ii0 + m\z)PR + m(z)Pi)e^-(-^*+P-"-)Vi;(Pr, z) 

= e^^(-^*+P-^-)[(7(70£; - ^tPt) + il^d, + m*{z)PR + m{z)PL\i>E{PT, z), 

(2.4) 

where a — +(— ) for positive(negative)-energy solution, pj- = {p^iP^)-, Xt — 
{x^,x'^), pt = \pt\ and 'JtPt = l^P^ + 7^P^- If we put E — E* cosh rj and 
Pt = E* sinh 77 with E* = ^^E'^ — p^, the Lorentz transformation of the spinor, 
represented by = e''''' ^^Z^?/', eliminates from (2.4) and from the Dirac equa- 
tion. Suppose that we do this Lorentz rotation for a fixed pj-. Then the Dirac 
equation is rewritten, in the second-order form, as 

[£;*2 + 52 _ |^(^)|2 ^ ,^/^(^)^3 _ rr^^{z)-i^^^\i>E{z) - 0, (2.5) 

where m(z) = mR{z) -|- imj{z), with both triR^z) and mj{z) being real-valued. 
Now let us introduce a set of dimensionless variables, using a parameter, a, which 
has a dimension of mass and whose inverse characterizes the thickness of the wall. 
Define 

■mR{z) = mof{az) = mo/(x), 

mi{z) mog{az) = mog{x), (2.6) 
X = az, T = at, e = E*/a, ^ = mo/a, 

where mo is the fermion mass in the broken phase. Then the Ansatz (2.4) and the 
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second-order equation (2.5) are expressed as 



'0(t, x) — e 



[(je7° + i7^— + ^f{x) - i^g{x)-f5]i^e{x), 



(2.7) 



and 



rf2 



+ - ^'(/(^)' + ^(^)') + i^f'i^h' - C9'{x)-f5^']Mx) = 0, (2.8) 



respectively. 

Our aim is to solve (2.8) and put it back into the Ansatz (2.7) to obtain 
the Dirac spinor, and to compute various currents from which the reflection and 
transmission coefficients are defined. As for f{x) and g{x), we do not specify their 
functional forms but only assume that 



and that \g{x)\ << 1, which means that the CP violation is small. (2.9) means 
that the system is in the broken phase at x ~ +00 while in the symmetric phase 
at x ~ —00. 

2.2. Outline of the DWBA to the Dirac equation 

Let us explain our strategy to construct the spinor solution. Since \g{x)\ is 
very small, we regard it as a perturbation and keep only quantities of 0{g^). First, 
put 




1, as X — > +00, 
0, as X — > —00, 






where ip^'^^x) is a solution to the unperturbed equation 




d2 





dx"^ 
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with an appropriate boundary condition. Then ijj^^^x) can be solved as 

^W{x)= j dx'G{x,x')V{x')ij^^\x'), (2.12) 

where 

V{x) = -ig\x)^^^\ (2.13) 
and G{x, x') is the Green's function satisfying 

[e' + ^ - enxf + iU'{x)^\^G^p{x, x') = -5^p5{x - x'), (2.14) 

with the same boundary condition as il}^^\x). Thus, to this order, the spinor 
solution to the Dirac equation is given by 

tlj{x) ~ e~*'^^^|[ae70+i73^+e/(a;)][^(°)(a;)+^(i)(a;)]-^e^(a;)75^(°Ha^)}- (2-15) 

2.3. Solution to the unperturbed equation 

In the following, any explicit forms of f{x) and g{x) are not necessary, except 
for the condition (2.9). We shall present general results on the wave function and 
currents below. The details of the calculation are given in Appendix. 

If we expand ■0^°^^) terms of the eigenspinors of 7^ as ■0^°^^) ~ <t>±{x)'^±i 
with 

7\^ = ±ra^ (s = l,2), (2.16) 



(t>±{x) must satisfy 



[6^ + ^ - ef{xf T U'{x)\<t>±{x) = 0. (2.17) 



Because of (2.9), the asymptotic forms of 0±(a;) should be 

(t>±{x) - 

where a — iyj — ^ and /3 = ie. Here we consider the case where the scattering 



e"^, e-"^ {x +00) 
e^^, {x —00) 



7 



states exist at x = +00; e > We denote the independent solutions to (2.17) as 



x) and 0^ = [(p±^'^\x)y which satisfy 



(2.18) 



(2.19) 



at X — > +00. Their asymptotic forms at x — > — 00 are 

(P^t'\x) ~ 7±(«, /9)e^" + 7±(«, -/3)e-^^^ 

</'i~"^(^) 7±(-«, /?)e^" + 7±(-«, -/9)e-^^ 

because of the symmetry of (2.17) under /3 1— > — Needless to say, ^±{01, P)* — 
7±(— a, — /3) and 7±(q;, — /3)* = 7±(— Z?)- With these solutions, the general solu- 
tion to (2.11) can be constructed as 



s=l,2 ± 



Multiplied by the operator, ae^^ + «7^^ + ^f{x), it becomes the general solution 
to the unperturbed (CP-conserving) Dirac equation, 

"fTe70 + ij^^ - Cf{x)]Mx) = 0. (2.20) 

As shown in Appendix, the coefficients, (Ai)*^"^"), are not needed to construct the 
general iPq{x). Eventually, the general solution to (2.20) is given by 

^Q{x)=[ae^'^ + zj'-^ + U-{x)]i;^'\x) 

= ^{4")M5r"^(x)i.i + (e + a)0L-'^)(x)<] (2.21) 

s 

+ Ai+\ae<f>^+''\x)ut + (^ - a)0L+"^(a;)<]}, 

where 

i^^'\x) = ^[4-Vi""^(^) + 4+V^"\^)]<- (2.22) 

s 

This ijj^^^x) is the starting point of the DWBA. Any boundary condition can be 
implemented by choosing and ^li"*""^ appropriately. 
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Now we impose a boundary condition to specify the scattering state. We 
consider a state in which the incident wave coming from a; = — oo (symmetric 
phase) is reflected in part at the bubble wall, while at x = +00 only the transmitted 
wave exists* This situation is achieved by setting = for cr = + and 

= for (7 = -, respectively, in (2.21) and (2.22). The Green's function 
(2.14) which matches this boundary condition can be explicitly constructed from 
the functions (f)^°'\x) and 0^ as shown in Appendix. Thus we are ready to 

compute the currents in the asymptotic regions. 



2.4. Wave function and currents in the asymptotic regions 

According to (2.15), the perturbed wave can be obtained once the unperturbed 
wave functions and the Green's function are known. We present only the results, 
while the derivations are given in Appendix. The asymptotic forms of the wave 
function are 



trans 



—e 



—ier+ax 



- h + -9 —00, 

eV 2q; 2^^ 



s 



{ip{x)a=+y 



—e 



-ieT+/3x 



e\2a 7+ (a, p) 2 



+ e l + (- - — —h--g{-oo) 

'- e V2q! 7+ (a, p) 2 



u_ 



g eV 2a 7+(a, -/?) 2 

^ \^ ^( \si( P 7-(-a,-/5) . 1 . 
+ + {-) -A- — ^TZ. - 7;9{-oo, 



e V 2a 7+(q;, — /?) 



(2.23) 



■k This boundary condition is referred to as Type II in ref. [8] . 
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and 



asym 



{a,f3,e,h ,l2)^{-a,-l3-eJl ,1^ 



(2.24) 



where 



/i= y dxg\x)(i>^I''\x)(i>^^''\x), 

—oo 
oo 

l2= j dxg\x)4)^^''\x)4)^^''\x). 



(2.25) 



From these, the asymptotic forms of the vector and axial- vector currents defined 

by 

3v = i^l^i^^ 3 a = V'T^'tsV', 

are (see Appendix) 

s 



OUn = 2e'M-,P)fj:U/^ ^[1 - a{-r'S 



where the corrections by the CP violation are 



2y?3^v 7+ («,/?) 



27^23^^ 7+(a,-/?) 



/2 + C.C.j, 
/2 + C.C.) . 



(2.26) 



(2.27) 
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The chiral currents defined by 



are 



A, 



A 



(^) 



(2.28) 



With these currents, one can compute the transmission and reflection coefli- 
cients, T^"^^ and R^'^\ respectively. From the vector currents, we have 



rp{a) 



trans 



(AT" 



1 + 



,.s+l 



A 



-aS' 



1 + 



Af^ 



Af^ 



where 



a 



/?|7+(«,/3)| 
7+(a, 



2 ' 



(2.29) 



(2.30) 



7+(a,/5) 

are those in the absence of the perturbation. As shown in Appendix, the unitarity 
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holds: 



T(-) + i?(-)=T(0) + i?W = l. 



(2.31) 



The differences in these coefficients between the particle and the anti-particle would 
vanish, even in the presence of CP violation, upon averaged over the initial mixed 
states such as the thermal equilibrium. This is because the difference has a factor. 



E.(- 



A 



A 



, in which A^g^ is representation-dependent and can take any value 



in a mixed state. Therefore, it is expected that there would be no net (vectorlike) 
particle current through the bubble wall. 

The transmission and reflection coefficients for the chiral fermion are defined 



as 



^L^L{R) ~ \JL{R),ct) \a^=0' yJ^cr) ' 

/ .3 \trans\ / f -3 \inc 

\3L{R),a) U?=o/ > 



(2.32) 



^L{R)^R{L) - \3R{L),a) I VL[R),u) ■ 

If we denote R{T) = i?(+)(T(+)) and R{f) = i?(-)(T(-)), we have 

= Tr^r 



Tl-^r 



Tr^r = Ti^L = 

Rl^r = Rr^l = 
Rr^l = Rl^r = 







2e|7+ 


>,/5)l' 






2e 7+ 


;«,/3)i' 






2e|7+ 


:«,/9)i' 







(2.33) 



2e|7+(a,/3)r 

7+(«, -0) 
7+(q;,/3) 

)'+' 



7+(q;,/3) 
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Among these, the following unitarity relations hold: 



Tl^l + Tl^r + Rl-*r = 1, 

(2.34) 

Tr^L + Tr^r + Rr-^l = 1- 

The difference between a particle and its anti-particle in the coefficients is, for 
example, 

Ai? = Rr^l - Rr^l = 2R^^\5''"' + S'^f^) = -2r(0)5*^^ (2.35) 

where we used 5'^'^'^ + R^^)^'^^^^ — (see Appendix). This suggests that we would 
have nonzero chiral-fermion current through the bubble wall even for the initial 
state in the equilibrium, if there is a nontrivial CP violation. (By 'nontrivial', we 
mean that g{x) does not proportional to f{x), since in that case the CP-angle can 
be removed by a constant bi-unitary transformation of the fermions, or explicitly, 
^tnc _ ^refl — as shown in Appendix.) 

3. Application to the Kink- Type Bubble Wall 

The prescription developed in the previous section is applicable to the situa- 
tion where the bubble formed in the first-order phase transition is moving with 
a constant velocity. The profile and motion of such a bubble would be deter- 
mined by analyzing the detail of dynamics of the phase transition. Far from the 
wall, deep in the symmetric/broken phase, the Higgs field is expected to have the 
vanishing/nonvanishing vacuum expectation value in the uniform vacuum. (This 
corresponds to the condition (2.9).) The simplest situation may be that in which 
the profile has the shape of the kink. This case, without CP violation, was studied 
in ref.[8]. Now we apply our procedure to that case, and evaluate the transmission 
and reflection coefficients. 
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Following [8], we suppose that the bubble- wall profile is given by 



1 + tanhx 

fix) = ^ ■ (3.1) 



Then the solutions to the Klein-Gordon-type equation, (2.17), are expressed in 
terms of the hypergeometric functions: 

^':%) = y-^'M - v f'Fi^^±^ + 1, -a + 1; ,), 

<i>^^^\y) = y-/\i - yf'Fi^^±l±^ + 1, ^±^,a + i-y) = (^^^ (,))*, 

(3.2) 

where we have changed the variable as y = (1 — tanhx)/2* These have the desired 
asymptotic behaviors as (2.18) and (2.19). Now 7±(q;,/3) are explicitly given by 

The transmission and reflection coefficients without the CP violation are 

^(0) ^ sm{na) sin(7r/3) 

sin[f(a + /5 + 0]sin[f(a + /3-0]' 
^ sin[f(a-/3 + 0]sin[f(«-/3-0] ^ " ^ 

sin[f(a + /3 + 0]sin[f(a + /3-0]' 



The effects of the CP violation can be evaluated, once the functional form of 
g{x) is given. One of the quantities we concern may be the difference in the reflec- 
tion/transmission coefficients between the particle and anti-particle. As mentioned 
above, the vectorlike one would vanish after averaging over the initial mixed state. 
However, the chiral one, which is given by (2.35), would survive even for the initial 
state in the equilibrium. 

★ The parameters a, (3 and ^ in ref.[8] are halves of ours. 
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We have executed numerical calculations for various g{x). Here we present the 
results for two typical cases in which g'{x) has different behaviors around the wall 
surface, i.e., g[x) = A6 ■ f{x)'^ and g{x) = A9-^f{x), where characterizes the 
magnitude of CP violation. We carried out the numerical calculation of AR at 
several energies. The a-dependence of AR/ AO for g{x) — A9- f{x)'^ is summarized 
in Table 1. The £'*-dependence of AR/A9 is plotted in Figs. 1 and 2 for the two 
choices of g{x) for several a. Prom these, we find the following: 

(1) The sign of AR varies not only with a but also with the functional form of 

(2) The absolute value of AR rapidly decreases as the energy of the incident 
particle increases. This feature may be independent of the choice of g{x), 
since the g'(x)-independent factor of AR is very small in the region far from 
E* ~ mo. 

(3) The degree of the decrease of \AR\ seems to be weakened as a increases, i.e. 
the bubble wall becomes thinner. 

4. Conclusions 

We have developed a general prescription to solve the Dirac equation in a 
bubble-wall background with CP violation, to the first order of CP violation. 
The transmission and reflection coefficients are obtained to this order, in terms 
of the functions of the energy, the mass in the broken phase and the wall thick- 
ness, which are extracted from the continuation of the solution to the second-order 
Klein-Gordon-like equation. It is found, irrespectively of the detailed form of the 
bubble-wall profile, that for the vector like current, the effect of CP violation will be 
canceled, when averaged over an initial mixed state such as thermal equilibrium, 
and that the net chiral current through the bubble wall would remain. Since the 
relation between the chemical potentials of the left- and right-handed quarks in 
the symmetric phase is different from that in the broken phase , this would leave 
net baryon number after the phase transition. 
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We have applied our prescription to the simple example in which the profile 
of the bubble wall is of the kink type. This situation may be realized when the 
phase transition proceeds slowly with degenerate minima of the free energy. Some 
numerical results are obtained in the case where the CP-violating angle has non- 
trivial spatial dependence. AR, which measures the difference in the reflection 
probabilities between the chiral fermion and anti-fermion, is evaluated for various 
energies and wall thickness, assuming some typical forms of g{x). The results sug- 
gest that irrespectively of the choice of g{x), |Ai?| rapidly decreases as the energy 
of incident fermions increases. Further the degree of decreasing seems to be larger 
as the wall becomes thicker. 

In order to determine whether and how much the net baryon number is left 
after the electroweak phase transition, we must know the detailed dynamics of 
the phase transition and particle distribution within both phases divided by the 
bubble wall. The baryon number left in the broken phase may be calculated by 
taking the average of the net fermion current through the bubble wall from the 
symmetric to the broken phase, and vice versa. Since the distribution function 
of the chiral fermions in the broken phase is different from that in the symmetric 
phase, as shown from the equilibrium analysis in [10], a quantity such as AR will 
play an important role in estimating the baryon asymmetry of the universe. As 
far as the cases we have studied are concerned, AR shows no drastic enhancement. 
On the other hand, even the sign of AR varies with the inverse wall-thickness a, as 
well as with the functional form of g{x). Thus whether the baryon or anti-baryon 
is left after the electroweak phase transition would be a subtle problem. We would 
hke to remark that |A^| does not need to be extremely small in so far as |gi(-|-oo)| 
is small enough to satisfy the experimental bound. In fact, in our second example 
{g{x) = A9f'{x)), g{+oo) vanishes irrespective of the value of A^. This fact might 
allow us to build a model which produces sufficient baryon number, by taking a 
suitable form of g{x). 

Our results will be applied to the scenario of the electroweak baryogenesis, 
which assumes the first-order phase transition. Although much should be revealed 



16 



about the mechanism of CP violation and the dynamics of the phase transition, 
such as the order of it and the motion of the bubble wall when it is of first order, 
the baryon asymmetry of the universe might be generated in the electroweak phase 
transition era if there is sufficient baryon number current through the bubble wall. 



APPENDIX 

In this appendix, we shall derive some equations used in the text. 
Unperturbed wave function 

The second-order unperturbed equation (2.17) is written as 

D^D±(l>i{x)=e^M^), (A.l) 

where is either or (j)^~'^\x) and the operator D± is defined by 

D± = ^^+Cf{x). (A.2) 

Multiplying (A.l) by D±, we have 

(D±D^)(D±0i+«)(x)) = 62(D±0(+")(a;)). 

Since the independent solutions to this equation are (l)^'^\x) and 0^ 
D±(f)^^'^\x) should be expressed as 

D±(P^+"\x) = c(+"V^"^(x) + c^'^cl^'^ix). 

Prom the behavior of the left-hand side at x ~ H-oo; 

D±<f>'^^"\x)r^{CTcy)e"^ 
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it follows that c. 



a and c. 



i-a) 



0. That is, 



In the Dirac representation, the 7-matrices are represented as 



(A.3) 



7 





-as 



7 



1 
-1 



and the eigenspinors of 7^ satisfying 



"^-^ -±i<, (s = 1,2) 
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(A.4) 



are expressed as 



u. 



1 

72 



( 1 \ 



\ / 



U4- = 



/ \ 

1 





(A.5) 



Then 7°ii^ = li^. The general solution to the unperturbed Dirac equation (2.20) 
is obtained by multiplying the most general ■0^°^(x) by the operator '^£7° + ^7^^ + 



+ (^;)(+")[ae0V+"^(a;)«i + - a)0L+"^(a;)M^] 
+ (^i)(+«)[(7e(/)L+"^(a;X + + a;)(/)^+")(x)iii] }. 



Here cre^^ '^\x)ut + {^ + a)(j)\_ "''{x)u^ and cre^Ii "''(a;)'u^ + (^ — q;)^!^ '^'{x)ut are 
not independent with each other, since the determinant of the coefficients of their 



(A.6) 
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components vanishes: 



— 0. Similarly the functions in the last two 



i — a ae 

lines of (A. 6) are not independent. Hence the most general solution is given by 
putting (741) = (74i)(+") = 0; the most general il)^'^\x) is 



V;(0)(a;) = + 4+"V^"^(a:)]«;. (A.7) 



(+a) ,(+a). 



The wave function satisfying the boundary condition which we adopted in the text 
is obtained by putting A^s = 0. 

Green's function 

Here we construct the Green's function satisfying 

[e^ + ^ - i'^f{xf + iif{x)^\^G^p{x, x') = l^a^G^p{x, x') = SapK^ - x'), 

(A.8) 

with the boundary condition discussed in the text. Put 



and introduce a unitary matrix defined by 



(A.9) 



Then, since A is diagonahzed as 

\ 

the Green's function is given by 

/ G+{x, x') 



A+ 



A-/ 



G{x,x) = U 



G-{x,x') 



G+{x,x') 



\ 



\ 



G-{x,x')J 
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where G±{x,x') satisfies 

A±G'±(x, x') = -5{x - x'), (A.IO) 

with the same boundary condition as G{x,x'). One can construct the Green's 
function G±{x, x') following the standard method'^'. If u±{x) {v±{x)) are solutions 
to A±u±(x) = 0(A±v±(x) = 0) and satisfy the boundary condition imposed on 
G{x, x') at X — > —oo{x — > +00), then the Green's function is expressed as 

{— ^ — \U±ix)v±ix'), ix < x') 

where W{u,v) is the Wronskian defined by 

W{u±,v±) — u±(x)-^v±(x) — v±{x)-^u±(x) = const. (^-12) 
ax ax 

For the boundary condition in question, there exists only right-moving wave at 
X ~ +00 so that we should choose, for each a, 

while we take the most general choice for u±{x); 

where is some constant. (The overall normalization is irrelevant, since the 
product of u and v is divided by the Wronskian in (A.ll).) By use of the asymptotic 
forms of (p±{x) at x ~ +00, 

W{u±, v±) — 2aa. 
Thus the Green's function G±{x,x') is 

I -U^^i'^'H^') + 4^V^'^"^(xO]0L+'^"^(a:), {x' <x) ' ' 
DWBA to the first order 
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As explained in Section 2.2, the function ipeix) in the Ansatz (2.7) is given by 

ipeix) = ilj^^\x) + i/j^^\x), (A.14) 

where ip^^^x) is a solution to the unperturbed equation with the boundary condi- 
tion above, that is, 

V'(°)(x)=^4'^V^""^(^)<, (A.15) 

s 

and ijj^^\x) is 

/O A 

Noting that in the Dirac representation, 75 ~ I ^ / ^^^^ 

and V{x) is 

the integrand of (A.16) is written as 

G{x,x')V{x')ilj^^\x') = -^G^^\x,x')g'{x')J2i-y^^f^(l'+''''\^'y'' 



u 



Therefore ^p^^\x) is 



s 



X 



|[0(_--«)(a;)+c(:V(_+""^x)] j dx'g'{x')<jt''''^{x')(l^''"^{x') 

X 

X 

+ 0L+""^(x) J dx'^'(x')[0L""''\x') + cL"Vi^"''^(^')]0^'"'\^')} 
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and its asymptotic forms are 



s 

X [(7-(-«, /5) + cL+^-(a, /3))e^-+(7-(-a, -/?) + cL+^-(a, -/3))e-^-] , 
for the positive-energy wave, and 



,s 



s Mi 



2q; 

s 

for the negative-energy wave, where 



— 00 

oo 

j dxg'{x)<p^^''\x)cjl^"\x). 



(A.17) 



With these asymptotic forms and those of ^p^^\x)^ 
(V'(°)(x))+°° = e^«-J]4"^<, 

s 

= 5]4")<[7+(<7a,/3)e/^- + 7+(<7a,-/3)e-/^-], 

s 

we can compute the asymptotic forms of the perturbed wave function using (2.15); 
^{x) = e-*^^^|[ae70+z73A+^/(^)][^(0)(^)+^(i)(^)]_,^^(^)^5^(0)(^)|_ (^.18) 
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For the positive-energy wave, we have 

(V'(x)^=+) '''''''' ^Q-i^r+ax ^ ^ 



(+) 
s 



X 



1 + 



+ e 



e 2a 



e 2q; 
h 



+ 



(+) 



X 



e2a 7+(q;, /3) 



X 



e 2q; 7+(a, — /?) 



e /3 7-(-«,-/3) j^ 

e 2q; 7_|_(a, — /9) 



u 



+ e 



1+ -r- 



U /3 7-(-«,-/9) 



e V 2a 7+(q!, — /3) 
The negative-energy wave is related to the positive one by 



bf^.-.(-oo))].i}. 

(A.19) 



When one of the relations (A. 3), 



asym 



{a,/3,e,h,l2)^{-a-p,-e,I^,i;) 



is evaluated at x ~ —00, it produces 



/37-(a, /3)e^^ - /?7-(«, -/Sje"^- = + a)[7+(a, /?)e^- + 7+(a, -/Jje" 



or 



/?7_(a,/3) = (C + a)7+(a,/3), 



(A.20) 
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These relations lead to 



2a 2a7+(a,/9) 2a7+(a,— /?) 

Hence all the factors including c^^^h in (A. 19) are the same, so that they are 
absorbed in the arbitrary coefficients A^^^ by redefining them as 



4« 



^ ^ + " (+) r ■ 
— ^ — c- h 
e 2a 



+ (-) 



e 2a 



1 + (-r|^i?(-oo) 



where the last factor of the left-hand side is put for later convenience. This gives 
the result (2.23). 

Currents and unitarity 

Now it is straightforward to calculate the asymptotic forms of the various currents, 
with the wave function (2.23). The straightforward calculation yields the same 
result as (2.26) except for (jy^^)^'''*"* and (j^ o-)*^**"*, which are 



■3 \ trans 



— —2iea 



EKf{ 



l + (-)V^ 



-g{+oo) + g{-oo) + (^^^^h + c.c.^ |, 



OX.r''' - E(-)1^^'^^ + [-9i+oo) + 9{-oo) + (^/i + C.C.)] }. 



Here one can show that 



2a 



Ji + c.c. = gi+oo) - g{-oo), 



(A.21) 



so that the result in (2.26) is obtained. Let us prove this equation. Define g{x) 
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g{x)—g{—(X)), theng'{x) = g\x) and5'(— oo) = 0. Integrating by parts, /i becomes 



J dx~g{x) 



From (A. 3), the first derivatives of 0±(a;) can be written as 



Hence, 



CXJ 

/i = ^(+00) - (e - a) j dx~g{x)[\cl>^^''\ 



X] 



which leads to 



2a 

+ a 



h + c.c. 



=^(+00). 



2a 



oc 

•)/ 



+ c.c. J / dxg{x 



Next we shall examine the unitarity. Without the perturbation, the unitarity 
relation. 



a 1 



7+(a, 



7+ 



(A.22) 



immediately follows from the current conservation which is derived from the Klein- 
Gordon- type equation, (2.17), or in other words, from the fact that the Wronskian 
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0+'"'') is independent of x. Evaluated at a; ~ +00 and at x ~ —00, the 
Wronskian is 



which is nothing but (A. 22). In the presence of the CP- violation, the sum of the 
transmission and reflection coefficients is, from (2.29), 



Ea(-) 



s+1 



4^) 



Now we can show that 



ginc ^ ^(0)^re/Z ^ 



(A.23) 



so that the unitarity of the perturbed wave holds. This relations is verified as 
follows. From the definitions. 



e 2a 



7-(-a,/3) 
7+(a,/3) 



7-(-",-/3) 



7+(a,-/5) 



+ c.c. 



Here the term in the parentheses is written as 



7_(-a,/5)^^(0)7-(-a,-/5) _ 7+(-q;, -/?)7-(-a, /3) + 7+(-a, /?)7-(-«, 



7+(a,/3) 



7+(a, -p) 



l7+(a,/?)r 



= 0, 



where we have used (A. 20) to rewrite 7_ in terms of 7_|_ in the denominator. 

The transmission and reflection coefficients of the chiral fermions are expressed 



as 



RR^L^R^^\l + S''"' + S''^f^), 
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from which it follows that 

Tl^l + Tl^r + Rl^r = 1 - (5'"' + R^'^'^S'^^') = 1, 
Tr^l + Tr^r + Rr^l = 1 + (5'"^ + = 1, 

because of (A. 23). 

When g(x) is a linear function of f{x), it is shown that 5*"''^ = 5"^^^^ — 0. This 
is verified as follows. By use of the differential equation (A.l), the integral in which 
g'{x) is replaced with f'{x) becomes 



72= / dxf'{x)(l)^^''\x)(l)^^''\x) 



/ dx 
2e J 

— oo 



dx'^ 



1 



/3 



[7+(a, /5)7-(a, - 7+(a, -/5)7_(a, 



Then, with the help of (A. 20), 



which is purely imaginary, so that (5'"^ = from the definition in (2.27), and 
^refl ^ because of (A. 23). 
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Table Caption 

Tablel. The values of AR/A9 for various energy of the incident particle and thickness 
of the bubble wall. We have chosen g{x) — A9f{x)^. The numerical values 
of E* and a are given in the unit of mo, the height of the wall. 

Figure Captions 

Figurel. AR/A6 as a function of E* for various a, in the case where g{x) — A6f{x)'^. 

The numerical values of E* and a are given in the unit of mo, the height of 
the wall. 

Figure2. AR/A9 as a function of E* for various a, in the case where g{x) = A9-^f{x). 

The numerical values of E* and a are given in the unit of mo, the height of 
the wall. 
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